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11.1 Introduction 

In this chapter we describe a selection of mathematical techniques and results that suggest inter- 
esting links between the theory of gratings and the theory of homogenization, including a brief 
introduction to the latter. By no means do we purport to imply that homogenization theory is 
an exclusive method for studying gratings, neither do we hope to be exhaustive in our choice 
of topics within the subject of homogenization. Our preferences here are motivated most of 
all by our own latest research, and by our outlook to the future interactions between these two 
subjects. We have also attempted, in what follows, to contrast the "classical" homogenization 



(Section 11.1.2), which is well suited for the description of composites as we have known them 



since their advent until about a decade ago, and the "non-standard" approaches, high-frequency 



homogenization (Section 11.2) and high-contrast homogenization (Section 11.3), which have 
been developing in close relation to the study of photonic crystals and metamaterials, which ex- 
hibit properties unseen in conventional composite media, such as negative refraction allowing 
for super-lensing through a flat heterogeneous lens, and cloaking, which considerably reduces 
the scattering by finite size objects (invisibility) in certain frequency range. These novel electro- 
magnetic paradigms have renewed the interest of physicists and applied mathematicians alike 
in the theory of gratings [HI. 



11.1.1 Historical survey on homogenization theory 



The development of theoretical physics and continuum mechanics in the second half of the 1 9th 
and first half of the 20th century has motivated the question of justifying the macrosopic view 
of physical phenomena (at the scales visible to the human eye) by "upscaling" the implied 
microscopic rules for particle interaction at the atomic level through the phenomena at the 
intermediate, "mesoscopic", level (from tenths to hundreds of microns). This ambition has 
led to an extensive worldwide programme of research, which is still far from being complete 
as of now. Trying to give a very crude, but more or less universally applicable, approximation 
of the aim of this extensive activity, one could say that it has to do with developing approaches 
to averaging out in some way material properties at one level with the aim of getting a less 
detailed, but almost equally precise, description of the material response. Almost every word 
in the last sentence needs to be clarified already, and this is essentially the point where one 
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could start giving an overview of the activities that took place during the years to follow the 
great physics advances of a century ago. Here we focus on the research that has been generally 
referred to as the theory of homogenization, starting from the early 1970s. Of course, even at 
that point it was not, strictly speaking, the beginning of the subject, but we will use this period 
as a kind of reference point in this survey. 

The question that a mathematician may pose in relation to the perceived concept of "aver- 
aging out" the detailed features of a heterogeneous structure in order to get a more homogeneous 
description of its behaviour is the following: suppose that we have the simplest possible linear 
elliptic partial differential equation (PDE) with periodic coefficients of period r\ > 0. What is 
the asymptotic behaviour of the solutions to this PDE as T] — )► 0? Can a boundary-value prob- 
lem be written that is satisfied by the leading term in the asymptotics, no matter what the data 
unrelated to material properties are? Several research groups became engaged in addressing 
this question about four decades ago, most notably those led by N. S. Bakhvalov, E. De Giorgi, 
J.-L. Lions, V. A. Marchenko, see [2l, [El, [ffl, iQ for some of the key contributions of that 
period. The work of these groups has immediately led to a number of different perspectives 
on the apparently basic question asked above, which in part was due to the different contexts 
that these research groups had had exposure to prior to dealing with the issue of averaging. 
Among these are the method of multiscale asymptotic expansions (also discussed later in this 
chapter), the ideas of compensated compactness (where the contribution by L. Tartar and F. 
Murat [l6l|, fT] has to be mentioned specifically), the variational method (also known as the 'T- 
convergence"). These approaches were subsequently applied to various contexts, both across 
a range of mathematical setups (minimisation problems, hyperbolic equations, problems with 
singular boundaries) and across a number of physical contexts (elasticity, electromagnetism, 
heat conduction). Some new approaches to homogenization appeared later on, too, such as the 
method of two-scale convergence by G. Nguetseng [8J and the periodic unfolding technique by 
D. Cioranescu, A. Damlamian and G. Griso [9|. Established textbooks that summarise these 
developments in different time periods, include, in addition to the already cited book IH, the 
monographs [jTOll . [[TT|. [[T2l. and more recently [fT3l . The area that is perhaps worth a separate 
mention is that of stochastic homogenization, where some pioneering contributions were made 
by S. M. Kozlov [14], G. C. Papanicolaou and S. R. S. Varadhan [15], and which has in recent 
years been approached with renewed interest. 

A specific area of interest within the subject of homogenization that has been rapidly de- 
veloping during the last decade or so is the study of the behaviour of "non-classical" periodic 
structures, which we understand here as those for which compactness of bounded-energy so- 
lution sequences fails to hold as T] — )■ 0. The related mathematical research has been strongly 
linked to, and indeed influenced by, the parallel development of the area of metamaterials and 
their application in physics, in particular for electromagnetic phenomena. Metamaterials can 
be roughly defined as those whose properties at the macroscale are affected by higher-order be- 
haviour as T] — 7- 0. For example, in classical homogenization for elliptic second-order PDE one 
requires the leading ("homogenised solution") and the first-order ("corrector") terms in the r\- 
power-series expansion of the solution in order to determine the macroscopic properties, which 
results in a limit of the same type as the original problem, where the solution flux ("stress" 
in elasticity, "induction" in electromagnetics, "current" in electric conductivity, "heat flux" in 
heat conduction) depends on the solution gradient only ("strain" in elasticity, "field" in elec- 
tromagnetics, "voltage" in electric conductivity, "temperature gradient" in heat condiction). If, 
however, one decides for some reason, or is forced by the specific problem setup, to include 
higher-order terms as well, they are likely to have to deal with an asymptotic limit of a different 



S. Guenneau et al.: Homogenization Techniques for Periodic Structures 11.3 

type for small 77 , which may, say, include second gradients of the solution in its constitutive law. 
One possible reason for the need to include such unusual effects is the non-uniform (in r\) ellip- 
ticity of the original problems or, using the language of materials science, the high-contrast in 
the material properties of the given periodic structure. Perhaps the earliest mathematical exam- 
ple of such degeneration is the so-called "double-porosity model", which was first considered 
by G. Allaire [HH and T. Arbogast, J. Douglas, U. Homung [[TTl in the early 1990s. A detailed 
analysis of the properties of double-porosity models, including their striking spectral behaviour 
did not appear until the work [ 18] by V. V. Zhikov. We discuss the double-porosity model and 



its properties in more detail in Section [TT3 



Before moving on to the next section, it is important to mention one line of research within 
the homogenization area that has had a significant role in terms of application of mathematical 
analysis to materials, namely the subject of periodic singular structures (or "multi-structures", 
see [ 19|). While this subject is clearly linked to the general analysis of differential operators on 
singular domains (see L20J), there has been a series of works that develop specifically homog- 
enization techniques for periodic structures of this kind (also referred to as "thin structures" in 
this context), e.g. [|2T]| . Il22l . It turns out that overall properties of such materials are similar 
to those of materials with high contrast. In the same vein, it is not difficult to see that com- 
pactness of bounded-energy sequences for problems on periodic thin structures does not hold 
(unless the sequence in question is suitably rescaled), which leads to the need for non-classical, 
higher-order, techniques in their analysis. 

11.1.2 Multiple scale method: Homogenization of microstructured fibers 

Let us consider a doubly periodic grating of pitch r\ and finite extent such as shown in Fig jll.l[ 
An interesting problem to look at is that of transverse electric (TE) modes — when the mag- 
netic field has the form (0,0,//) — propagating within a micro-structured fiber with infinite 
conducting walls. Such an eigenvalue problem is known to have a discrete spectrum: we look 
for eigenfrequencies (O and associated eigenfields H such that: 



'^): 



^^^"^"'-^ ""<'"/■ 

where we use the convention x = {x\,X2), dQ.f denotes the boundary Q.f, and n = {n\,n2) is 
the normal to the boundary. Here, £o/io = c^^ where c is the speed of light in vacuum and we 
assume that matrix coefficients of relative permittivity £;j(y), with ?, j = 1 , 2, are real, symmetric 
(with the convention y == {y\-,y2)), of period 1 (in yi eiy2) and satisfy: 

M| ||2 > e,j{y)^,^j > m| ||2 , V| G IR2 , Vy G y = [0, 1]^ , (11.1) 

where | ^ | = {t,^ + ^,2), for given strictly positive constants M and m. This condition is met 
for all conventional dielectric medi4B 



'When the periodic medium is assumed to be isotropic, £ij{y) — e(y)5,y, with the Kronecker symbol 5,; = 1 
if / = j and otherwise. For instance, ( |ll.l| i has typically the bounds M = 13 and m — 1 in optics. One class of 
problems where this condition ( |ll.l| i is violated (the bound below, to be more precise) is considered in Section 
|11.3|on high-contrast homogenization. 
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Y2=X2/r| 




Figure 11.1: A diagram of the homogenization process: when the parameter rj gets smaller (rj < rj'), the number 
of cells inside the fixed domain Q,f becomes larger When 7] ^ 1, Q.f is filled with a large number of small 
cells, and can thus be considered as an effective (or homogenized) medium. Such a medium is usually described 
by anisotropic parameters depending upon the resolution of auxiliary ("unit cell") problems set on the reseated 
microcopic cell Y which typically contains one inclusion D. 



We can recast (^n) as follows: 






dxi 



with 



a'Wx))=^-fjV"« 



^'J \riJ dxj ' 
The multiscale method relies upon the following ansatz: 

H = //o(x) + 77//i(x,y) + 7]2//2(x,y) + ... 



(11.2) 



where ///(x, y), / = 1, 2, ... is a periodic function of period Y in y. 

In order to proceed with the asymptotic algorithm, one needs to rescale the differential operator 

as follows 



+ ■ 



dxi V (^Zi dji ) V ^^i ^yi ) 

where d jdzi stands for the partial derivative with respect to the /th component of the macro- 
scopic variable x. 
It is useful to set 
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a'{H) = a^ + 7]a| + ri^o^ + ... 

what makes ( |11.3p more compact. 

Collecting coefficients sitting in front of the same powers of r/ , we obtain: 

and so forth, all terms being periodic in y of period 1. 
Upon inspection of problem (^tj), we gather that 



so that at order r\ 



and at order ry 



-1 



K):-|-aj = 






(the equations corresponding to higher orders in T] will not be used here). 

Let us show that (J^) provides us with an equation (known as the homogenized equation) 

associated with the macroscopic behaviour of the microstructured fiber. Its coefficients will be 

obtained thanks to {£/) which is an auxiliary problem related to the microscopic scale. We will 

therefore be able to compute Hq and Hi thus, in particular, the first terms of// and a'. 

In order to do so, let us introduce the mean on Y, which we denote < . >, which is an operator 

acting on the function g of the variable y: 



< 



^>= ]Y]J J Siyuy2)dyidy2 



where 1 1' | is the area of the cell Y. 

Applying the mean to both sides of (J^), we obtain: 



d : d ■ (0^ 

dzi dyi c^ 



< (^) >: -^^ < C7^ > - < 3^a| >= -^Hq < 1 > 



where we have used the fact that < . > commutes with d /dzi- 
Moreover, invoking the divergence theorem, we observe that 

where n = (ni , ^2) is the unit outside normal to dY of Y . This normal takes opposite values on 
opposite sides of 7, hence the integral over dY vanishes. 
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Altogether, we obtain: 



2 



<(^)>:-^<a^>=^//o 



dz 



c 



which only involves the macroscopic variable x and partial derivatives d /dzt with respect to the 
macroscopic variable. We now want to find a relation between < Oq> and the gradient in x of 
Hq. Indeed, we have seen that 






which from (s/) leads to 



We can look at (^1) as an equation for the unknown i/i(x,y), periodic of period ^ in y and 

parametrized by x. Such an equation is solved up to an additive constant. In addition to that, 

the parameter x is only involved via the factor dHo/dzj. Hence, by linearity, we can write the 

solution Hi (x, y) as follows: 

dHo{x) ■ 
^i(x,y =^ vv' y , 

where the two functions w-^(y), j = 1,2 are solutions to (^1) corresponding to dHo/dzj{x), 
j = 1,2 equal to unity with the other ones being zero, that is solutions to: 



^-^y--U^r^^^'^)-Mifi''^^^) 



with >v^(y), k = 1,2 periodic functions in y of period Y R 
Since the functions vv^(y) are known, we note that 



dHo ^ dHo <9w*(y) 



' ^ '^^ '^ ^^^ [ dzj + dy, J ''J ^^> [ dzj ' dz, dyj 
which can be written as 



cx^(x,y)^(%ny)+.-(y)^) 



-U ,dw''{y)\ dHo{x) 



dyj J dzk 
Lets us now apply the mean to both sides of this equation. We obtain: 

;• ^ N _ -1 dHo{x) 
< CTq > (Xj - Chom,*"^^ ' 

which can be recast as the following homogenized problem: 



^We note that {s/2) are two equations which merely depend upon £: ■ ' (y), that is on the microscopic properties 
of the periodic medium. The two functions w^ (defined up to an additive constant) can be computed once for all. 



independently of Q.f. 
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Figure 11.2: Potentials Vx (left) and Vy (right): The unit cell contains an elliptic inclusion of relative permittivity 
(£ = 4.0 + 3i) with minor and major axis a ~ 0.3 and b = 0.4 in silica (e = 1 .25). 



i' 



,-1 



dHQ{x) 



dxi V ''°'"''^ dx. 






,on dClf , 



.-1 



where e^^^ ^.^ denote the coefficients of the homogenized matrix of permittivity given by: 



.-1 



^/.((^^^■'^'-^^■'^'^)''- 



(11.4) 



As an illustrative example for this homogenized problem, we consider a microstructured waveg- 
uide consisting of a medium with relative permittivity e = 1 .25 with elliptic inclusions (of minor 
and major axes 0.3 cm and 0.4 cm respectively) with center to center spacing d = 0.1cm with 
an infinite conducting boundary i.e. Neumann boundary conditions in the TE polarization. 

We use the COMSOL MULTIPHYSICS finite element package to solve the annex prob- 



lem and we find that Chom from (11.4) writes as [26] 



1.9296204 -1.053308310-16 

-44.41744410-1^ 2.1127643 

with < e >Y= 2.2867255. The off diagonal terms can be neglected. 

If we assume that the transverse propagating modes in the metallic waveguide have a small 
propagation constant 7^1, the above mathematical model describes accurately the physics. 



We show in Fig. 1 1.3 a comparison between two TE modes of the microstructured waveguide 
and its associated anisotropic homogenized counterpart. Both eigenfrequencies and eigenfields 
match well (note that we use the waveguide terminology wavenumber k = ^J (iP- jc^ — y^). 



11.1.3 The case of one-dimensional gratings: Application to invisibility cloaks 

There is a case of particular importance for applications in grating theory: that of a periodic 
multilayered structure. Let us assume that the permittivity of this medium is £ = a in white 



layers and /3 in yellow layers, as shown in Fig. 1 1 .4 
Equation (£/2) takes the form: 



(^3): 



d 
dy 



-'<-<) 



(^-■'^' 
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Figure 11.3: Comparison between transverse electric fields TE21 and TE^i of a microstructured metallic waveg- 
uide for a propagation constant 7 = 0. Icm^ (wavenumbers k = 0.1101 cm^ and k = 0.5478cot^^ respectively), 
see left panel, with the TE21 and TEj,i modes of the corresponding homogenized anisotropic metallic waveguide 
forY=0.\cm^^ (k = 0.7607cm"' andk = 0.5201cm"^ where k = ^Op- jc'^ — "f- = ^OJ^Co/Xo — 7^ were obtained 
from the computation of eigenvalues CO of homogenized problem {:3^q)), see right panel 



with w{y) , periodic function in y of period 1 . 

We deduce that 

dw 



Noting that 



dw 
Y dy 



w(l)-w(0) 



dy 
0, this leads to 



\+Ce{y) 



J^{l+C£{y))dy = 0. 



Since I 7 1= 1, we conclude that 



C=-<£> 

The homogenized permittivity takes the form: 

^-1 



^/(.-w..-w^)..v 



"'hom I Y 

= <£-\y)> - <:£-\y)+C> 

=< £^\y) > - <£^\y)> + << £{y) >^ >= < £{y) > 



-1 



We note that if we now consider the full operator i.e. we include partial derivatives in yi and y2, 
the anisotropic homogenized permittivity takes the form: 

<£{y)-' 



> 







^hom I <£{y)>'^ 
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Figure 11.4: Schematic of iiomogenization process for a one- dimensional grating with homogeneous dielectric 
layers of permittivity a and /3 in white and yellow regions. When r\ tends to zero the number of layers tends to 
infinity, and their thicknesses vanish, in such a way that the width of the overall stack remains constant. 



as the only contribution for e^^^ j^ is 1/ | 7 | /y e^^(y) Jy. 

As an illustrative example of what artificial anisotropy can achieve, we propose the design of an 
invisibility cloak. For this, let us assume that we have a multilayered grating with periodicity 
along the radial axis. In the coordinate system (r, 0), the homogenized permittivity clearly has 
the same form as above. If we want to design an invisibility cloak with an alternation of two 
homogeneous isotropic layers of thicknesses J^ and dg and permittivities a, /3, we then need to 
use the formula 

er 1 + 77 V« PA 1 + ^ ' 

where r\ = dg/dA is the ratio of thicknesses for layers A and B and J4 + Jr = 1. 



We now note that the coordinate transformation / = i?i + r ^^ ' can compress a disc 
r < R2 into a shell Ri < r < R2, provided that the shell is described by the following anisotropic 
heterogeneous permittivity [|27ll e'^'°* (written in its diagonal basis): 



£r 



cloak 



Ri 



r' -Ri 



., cloak 
-6 



Ri 



(11.5) 



R2-R\J \ r' J ' " \R2-R1 

where Ri and R2 are the interior and the exterior radii of the cloak. Such a metamaterial can be 



approximated using the formula ( 1 1.5 1, as first proposed in [i28il . which leads to the multilayered 



cloak shown in Fig. 1 1 .5 



11.2 High-frequency homogenization 

Many of the features of interest in photonic crystals [|44ll45]| . or other periodic structures, such as 
all-angle negative refraction fi46ll47ll48ll49ll or ultrarefraction [l50l[5T]| occur at high frequencies 
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Figure 11.5: Propagation of a plane wave of wavelength 7 lO^^m (red in the visible spec- 
trum) from the left on a multilayered cloak of inner radius Ri = 1.5 10^ m and outer ra- 
dius R2 = 3 10^ m, consisting of 20 homogeneous layers of equal thickness and of respective 
relative permittivities 1680.70,0.25,80.75,0.25,29.39,0.25,16.37,0.25,10.99,0.25,8.18,0.25,6.50,0.25,5.40, 
0.25,4.63,0.25,4.06,0.25 in vacuum. Importantly, one layer in two has the same permittivity. 



where the wavelength and microstructure dimension are of similar orders. Therefore the con- 
ventional low-frequency classical homogenisation clearly fails to capture the essential physics 
and a different approach to distill the physics into an effective model is required. Fortunately 
a high frequency homogenisation (HFH) theory as developed in [ITTl is capable of capturing 
features such as AANR and ultra-refraction [52] for some model structures. Somewhat tan- 
gentially, there is an existing literature in the analysis community on Bloch homogenisation 
[|53ll54ll55l[56l . that is related to what we call high frequency homogenisation. There is also 
a flourishing literature on developing homogenised elastic media, with frequency dependent 
effective parameters, based upon periodic media [|38l . There is therefore considerable inter- 
est in creating effective continuum models of microstructured media that break free from the 
conventional low frequency homogenisation limitations. 



11.2.1 High Frequency Homogenization for Scalar Waves 

Waves propagating through photonic crystals and metamaterials have proven to show different 
effects depending on their frequency. The homogenization of a periodic material is not unique. 
The effective properties of a periodic medium change depending on the vibration modes within 
its cells. The dispersion diagram structure can be considered to be the identity of such a mate- 
rial and provides the most important information regarding group velocities, band-gaps of dis- 
allowed propagation frequency bands, Dirac cones and many other interesting effects. The goal 
of a homogenization theory is to provide an effective homogeneous medium that is equivalent, 
in the long scale, to the initial non-homogeneous medium composed of a short-scale periodic, 
or other microscale, structure. This was achieved initially using the classical theory of homog- 
enization H |34l [m |35l [SGl and yields an intuitively obvious result that the effective medium's 
properties consist of simple averages of the original medium's properties. This is valid so long 
as the wavelength is very large compared to the size of the cells (here we focus on periodic 
media created by repeating cells). For shorter wavelengths of the order of a cell's length a more 
general theory has been developed [37J that also recovers the results of the classical homoge- 
nization theory. For clarity we present high frequency homogeniaztion (HFH) by means of an 
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illustrative example and consider a two-dimensional lattice geometry for TE or TM polarised 
electromagnetic waves. With harmonic time dependence, exp(— /^?) (assumed understood and 
henceforth suppressed), the governing equation is the scalar Helmholtz equation, 

V\ + Q}u = 0, (11.6) 

where u represent Ez and Hz, for TM and TE polarised electromagnetic waves respectively, and 
Q^ = rp-ap- /c^. In our example the cells are square and each square cell of length 2/ contains 
a circular hole and the filled part of the cell has constant non-dimensionalized properties. The 
boundary conditions on the hole's surface, namely the boundary ^^2, depend on the polarisa- 
tion and are taken to be either of Dirichlet or Neumann type. This approach assumes infinite 
conducting boundaries which is a good approximation for micro- waves. We adopt a multiscale 
approach where I is the small length scale and L is a large length scale and we set r\ = l/L <^ 1 
to be the ratio of these scales. The two length scales let us introduce the following two indepen- 
dent spatial variables, ^t = Xj/l and Xt = xi/L. The cell's reference coordinate system is then 



— 1 < (^ < 1 . By introducing the new variables in equation ( 1 1 .6 1 we obtain, 

M(X,|),^,.^,+f22M(X,|)+27]M(X,^),^,.;,,.+7]2M(X,|),x,X, = 0. (11.7) 

We now pose an ansatz for the field and the frequency. 



i<(X,^)=i/0(X,.§) + T]l<i(X,|)+T]2i,2(X, .§) + ..., 

Q^ = al + r\Q\ + r\'^al + ... (11.8) 

In this expansion we set ^o to be the frequency of standing waves that occur in the perfectly pe- 
riodic setting. By substituting equations (11.8) into equation ( |11.7[ ) and grouping equal powers 



of e through to second order, we obtain a hierarchy of three ordered equations: 

"o,§,^,+^o"o = 0, (11.9) 

"i,§,-.^, + ^o"i = -2mo,^,x, - ^i"o, (11.10) 

"2,<^,-^, + ^0"2 = -M0,X,X, -2Mi^^.x, -^l"l -^2"0- (11-11) 

These equations are solved as in [|40l l37ll and hence the description is brief. 

The asymptotic expansions are taken about the standing wave frequencies that occur at 



the corners of the irreducible Brillouin zone depicted in Fig. 1 1.6 It should be noted that not 
all structured cells will have the usual symmetries of a square, as in Fig. |11.6| ^a) where there 
is no reflexion symmetry from the diagonals. As a consequence the usual triangular region 
YXM does not always represent the irreducible Brillouin zone and the square region TMXN 
should be used instead. Also paths that cross the irreducible Brillouin zone have proven to yield 
interesting effects namely along the path MX' for large circular holes [39|. 

The subsequent asymptotic development considers small perturbations about the points 
r, X and M so that the boundary conditions of u on the outer boundaries of the cell, namely 
dSi, read, 

"l^;=i = ±" 1^,-1 and M^.|^_i =±M^^.|^__i, (11.12) 

where the -I-, — stand for periodic and anti-periodic conditions respectively: the standing waves 
occur when these conditions are met. The conditions on dS2 are either of Dirichlet or Neumann 
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Figure 11.6: Panel (a) An infinite square array of split ring resonators with the elementary cell shown as the 
dashed line inner square. Panel (b) shows the irreducible Brillouin zone, in wavenumber space, used for square 
arrays in perfectly periodic media based around the elementary cell shown of length 21 (I — I in (b)). Figure 
reproduced from Proceedings of the Royal Society f40^]. 



type. The theory that follows is similar for both boundary condition cases, but the latter one 
is illustrated herein. Neumann boudary condition on the hole's surface or equivalently electro- 
magnetic waves in TE polarization yield, 



du I 

which in terms of the two-scales and a,(X, ^ ) become 



0. 



Uo,^i'^i = 0, {Uofo^xi + ui^^.)ni = 0, (ui^xi + U2,^i)ni = 0. 



(11.13) 



(11.14) 



The solution of the leading order equation is by introducing the following separation of variables 
"0 = /o(X)t/o(^;^o)- It is obvious that /o(X), which represents the behaviour of the solution 
in the long scale, is not set by the leading order equation and the resulting eigenvalue problem is 
solved on the short- scale for Qq and Uq representing the standing wave frequencies and the as- 
sociated cell's vibration modes respectively. To solve the first order equation ( |11.10[ ) we take the 
integral over the cell of the product of equation (|11.10|) with Uq minus the product of equation 



( 1 1 .9 ) with Ui/fo and this yields Qi = 0. It then follows to solve for mi (X, ^ ) = /o,;c,. (X) t/i,. ((§ ] 



where the vector Ui is found as in [I40II . By invoking a similar solvability condition for the 
second order equation we obtain a second order PDE for /o(X), 



Tijfo,XiXj + ^2/0 = where, 



'•J 



IfsU^dS 



for 



I, J 



1,2 



(11.15) 



entirely on the long scale with the coefficients Tjj containing all the information of the cell's 
dynamical response and the tensor tij represents dynamical averages of the properties of the 
medium. For Neumann boundary conditions on dS2 its formulation reads, 

tu = j juldS + j j{U,^^^Vo-Ui^UQ^^;)dS for /=lor2, (11.16) 
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Figure 11.7: Tire dispersion diagram for a doubly periodic array of square cells with circu lar inclusions, of radius 
0.4, free at their inner boundaries shown for the irreducible Brillouin zone of Fig. \11.6\ The dispersion curves 
are shown in solid lines and the asymptotic solutions from the high frequency homogenization theory are shown in 
dashed lines. Figure reproduced from Proceedings of the Royal Society l\4d\l . 



([/i.,^,.C/o-f/i^.C/o,^,.)J5 for i^j. 



(11.17) 



Note that there is no summation over repeated indexes for tu. The tensor depends on the bound- 
ary conditions of the holes and has a different form if Dirichlet type conditions are applied on 

^52. 

The PDE for /o has several uses, and can be verified by re-creating asymptotically the 
dispersion curves for a perfect lattice system. One important result of equation ( 11.15[ ) is its 



use in the expansion of Q. namely in equation (11.8). In order to obtain Q.2 as a function of 
the Bloch wavenumbers we use the Bloch boundary conditions on the cell to solve for /o(X) = 
exp{iKjXj/ri), where Kj = Kj — dj with dj = 0, ;r/2, — ;r/2 depending on the location in the 
Brillouin zone. The asymptotic dispersion relation now reads, 



^~^o + 



la. 



-KiKj. 



(11.18) 



Equation (11.18) yields the behaviour of the dispersion curves asymptotically around the stand- 
ing wave frequencies that are naturally located at the edge points of the Brillouin zone. Fig. 



11.8 illustrates the asymptotic dispersion curves for the first six dispersion bands of a square 



cell geometry with circular holes. 
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An assumption in the development of equation ( 1 1 . 1 8 1 is that the standing wave frequen 



cies are isolated. But one can clearly see in Fig. 1 1 .7 that this is not the case for third standing 
wave frequency at point T as well as for the second standing wave frequency at point X. A small 
alteration to the theory [|40ll enables the computation of the dispersion curves at such points by 
setting, 

UQ = f^\X)ui;\^-^Q) (11.19) 

where we sum over the repeated superscripts (/). Proceeding as before, we multiply equation 



(ll.lOlby Uq , %\xhs\xdici u\{{U^J)^. + Q^Uq) then integrate over the cell to obtain. 



(-l_Xj^i + a\^^^f^^=Q, for m-l,2,...,p (11.20) 

Qi is not necessarily zero, and 

P^jmi = ///f/,S"^<i -<J.f/,S'V5, B„, = JJu^'X""^dS. (11.21) 

There is now a system of coupled partial differential equations for the /q ' and, provided 
Q.I 7^ 0, the leading order behaviour of the dispersion curves near the Q.q is now linear (these 
then form Dirac cones). 

For the perfect lattice, we set /g = /q exp{iKjXj/ri) and obtain the following index 
equations, 

{i'^Aj„,, + a\B„,i)f^'^=0, for m=l,2,...,p (11.22) 



The system of equation ( |11.22| ) can be written simply as, 

CFo = 0, (11.23) 

with C// = ^iB// and C,„/ = iKjAj,„i/r\ for I ^ m. One must then solve for Q}^ = ±y/aijlqK]/ri 
when the determinant of C vanishes and insert the result in, 

1 



^^^o±:^7^^/aijKiKj. (11.24) 

ZhiQ 

If the ^1 are zero one must go to the next order. 

11.2.1.1 Repeated eigenvalues: quadratic asymptotics 

If Q.I is zero, ui = /qx ^i (we again sum over all repeated (/) superscripts) and we advance 



to second order using ( |11.11[ ). Taking the difference between the product of equation ( |11.11| ) 
with Uq and M2(f^o.£,<?, + ^o^o) and then integrating over the elementary cell gives 

+a2yWjj^f;Hf;Wj5^0, for m=l,2,...,;? (11.25) 

as a system of coupled PDEs. The above equation is presented more neatly as 

fojCiXi^mi + fo%xPkjmi + ^ifo^mi = 0, for m = 1 , 2, . . . , p. (11 .26) 
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Figure 11.8: THe dispersion diagrams for ci doiihly periodic array' of sq uare c ells wilKsplit ring inclusions, free at 
their inner boundaries shown for the irreducible Brillouin zone of Fig. \11.6\ The dispersion curves are shown in 
solid lines and the asymptotic solutions from the high frequency homogenization theory are shown in dashed lines. 
Figure reproduced from Proceedings of the Royal Society 1 40] . 



For the Bloch wave setting, using /q ^(X) = /q ' exp{iKjXj/ri) we obtain the following system, 



and this determines the asymptotic dispersion curves. 



^D^,w + a^B,„, US'^=0, for m=l,2,...,/7 



(11.27) 



11.2.1.2 The classical long wave zero frequency limit 

The current theory simplifies if one enters the classical long wave, low frequency limit where 
Q? ~ O(e^) as Uq becomes uniform, and without loss of generality is set to be unity, over the 



elementary cell. The final equation is again (11.15) where the tensor tij simplifies to 

tu = j j^dS + j JU.^^^AS, Uj = j JU,.^^AS for i^j (11.28) 

(with no summation over repeated suffices in this equation) and Tij = tij/ J JgdS. 



11.2.2 Illustrations for Tranverse Electric Polarized Waves 



Let us now turn to some illustrative examples. We present in Fig. 11.8 the TE polarization 
waves for three types of SRR's (Split Ring Resonator's). Equation ( |11.15[ ) represents the wave 
propagation in the effective medium. It is noticable that the Tjj coefficients depend on the stand- 
ing wave frequency and that Tn is not necessarily equal to T22 in order to yield an anisotropic 
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effective medium for each separate frequency. Near some of the standing wave frequencies the 
anisotropy effects are very pronounced and well explained by the no longer elliptic equation 



(11.15). 



In the above equations U\. is a solution of, 



Ui,M> = ^^ (11-29) 

with boundary conditions (/o.x, + "i,£,)^! = on the hole boundary. If the medium is homoge- 



neous as it is in the illustrative examples herein, equation ( 1 1 .29 1 is the same as that for Uq, but 
with different boundary conditions. The specific boundary conditions for Ui . are 



U 



1 ■ £•"/ 



n; 



for J = 1,2, 



(11.30) 



where n, represent the normal vector components to the hole's surface. The role of Ui is to 
ensure Neumann boundary conditions hold and the tensor contains simple averages of inverse 
permittivity and permeability supplemented by the correction term which takes into account the 
boundary conditions at dS2- Equation (11.28) is the classical expression for the homogenised 



coefficient in a scalar wave equation with constant material properties; (11.29) is the well- 
known annex problem of electrostatic type set on a periodic cell, see [l4l[TT]|, and also holds 
for the homogenised vector Maxwell's system, where Ui now has three components and /, j = 

1,2,3 igB mil. 



11.2.2.1 Cloaking in metamaterials 



SRRs with 4 holes are now used and the dispersion diagrams are in Fig. 1 1.8 (b). The flat band 
along the MT path is interesting for the fifth mode and we choose to illustrate cloaking effects 
that occur here. In Fig. 11.9;a), we set an harmonic source at the corresponding frequency 
Q. = 2.8 in an 8 X 8 array of SRRs and observe a wave pattern of concentric spherical modes. 
As can be seen in Figs. |11.9|^b) and|11.9|^c) a plane wave propagating at frequency Q = 2.8 



demonstrates perfect transmission through a slab composed of 38 SRRs but also cloaking of a 
rectangular inclusion where no scattering is seen before or after the metamaterial slab. Panel (d) 



of Fig. 1 1 .9 shows the location in the band structure that is responsible for this effect. Note that 
the frequency of excitation is just below the Dirac cone point located at Q = 2.835 where the 
group velocity is negative but also constant near that location of the Brillouin zone illustrated 



through an isofrequency plot of lower mode of the Dirac point in Fig. 1 1 .9 e). In constrast with 
the isotropic features of panel (e), those of panels (f) and (g) show ultra-flattened isofrequency 
contours that relate to ultra-refraction, a regime more prone to omni-directivity than cloaking. 



The asymptotic system of equations ( 1 1 .20 ) describing the effective medium at the Dirac point 
can be uncoupled to yield one same equation for all /q 's. 










(11.31) 



After some further analysis, the PDF for /q is responsible for the effects at the frequency 
chosen fl = 2.8. 



11.2.2.2 Lensing via AANR and St Andrew's cross in metamaterials 

We observe all- angle-negative-refraction effect in metamaterials with SRRs with 8 holes. The 



dispersion curves in Fig. 11.8 c) are interesting, as the second curve displays the hallmark of 
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Figure 11.9: Cloalcing in square arrays of SRRs witli four Iwles: A source at fre quency Q. — 2.8, located in 



tlie center of a square metamaterial consisting of 64 SRRs sliaped as in Fig. 11.8\ b) produces a wave pattern 
reminiscent of (a) concentric spherical field, (b) cloaking of a rectangular inclusion inside a slab of a metamaterial 
consisting of 38 SRRs and (c) scatte ring o f a plane wave from the same rectangular hole as the previous panel, (d) 
Zoom in dispersion diagram of Fig. \11.8\ b). Panels (e), (f) and (g) present isofrequency plots of the respective the 
lower, middle and upper modes of the Dirac point. Figure reproduced from Proceedings of the Royal Society i l4Ql/ . 
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(b) 





Figure 11.10: Lensing via AANR and St Andrew's cross in square arrays of SRRs with eiglit lioles: (a) A line 
sourc e at frequency D, — 1.1375 located above a rectangular metamaterial consisting of of 90 SRRs as in Fig. 
[TXSJ'cj displays an image underneath (lensin g); (b ) A line source at frequency D. = 1.25 located inside a square 



metamaterial consisting of 49 SRRs as in Fig. \ll.W{c) displays the dynamically induced anisotropy of the effective 
medium; (c) Zoom in dispersion diagram of Fig, \ll.8\ c). Note that each cell in the arrays in (a) and (b) has been 
rotated through an angle ;r/4. Figure reproduced from Proceedings of the Royal Society l\40il . 



an optical band for a photonic crystal (it has a negative group velocity around the F point). 
However, this band is the upper edge of a low frequency stop band induced by the resonance 
of a SRR, whereas the optical band of a PC results from multiple scattering, which thus arises 
at higher frequencies. We are therefore in presence of a periodic structure behaving somewhat 
as a composite intermediate between a metamaterial and a photonic crystal. One of the most 
topical subjects in photonics is the so-called all-angle-negative- refraction (AANR), which was 
first described in P6ll . AANR allows one to focus light emitted by a point, onto an image, 
even through a flat lens, provided that certain conditions for AANR are met, such as convex 
isofrequency contours shrinking with frequency about a point in the Brillouin zone [49]. In Fig. 
11.10[ we show such an effect for a perfectly conducting photonic crystal (PC) in Fig. 11.10 a). 
In order to achieve AANR, we choose a frequency on the first dispersion curve (acoustic band) 
in Fig. 11.8 c), and we take its intersection with the light line Q =| K | along the XT path. This 



means that we achieve negative group velocity for waves propagating along the XT direction 
of the array, hence the rotation by an angle tt/A of every cell within the PC in panel (b) of Fig. 
11.10[ This is a standard trick in optics that has the effect of moving the origin of the light-line 
dispersion to X as, relative to the PC, the Bloch wavenumber is along XT. This then creates 
optical effects due to the interaction of the light-line with the acoustic branch, this would be 
absent if F were the light-line origin. 



The anisotropy of the effective material is reflected from coefficients Tn = —5.53 and 
T22 = 0.2946. The same frequency of the first band is reachable at point A^ of the Brillouin 
zone. By symmetry of the crystal, we would have Tn = 0.2946 and T22 = —5.53. The resultant 
propagating waves would come from the superposition of the two effective media described 



above. Fig. ll.lO^b) illustrates this anisotropy as the source wave only propagates at the pre- 



scribed directions. 
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Figure 11.11: For the two dimensional example we sliow tiie geometry of tlie doubly periodic simply supported 
plate (the dots represent the simple supports) in panel (a) with the elementary cell shown by the dotted lines and 
in (b) the irreducible Brillouin zo ne w ith the lettering for wavenumber positions shown. Figure reproduced from 
Proceedings of the Royal Society 4591/ . 
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Figure 11.12: The dispersion diagra m for a doubly periodic array of point simple supports shown for the irre- 
ducible Brillouin zone of Fig. 11.11 The figure shows the dispersion curves as solid lines. As dashed lines, the 
asymptotic solutions from the highfrequency homogenization theory are shown. Figure reproduced from Proceed- 
ings of the Royal Society 4591/ 



11.2.3 Kirchoff iMve Plates 

HFH is by no means limited to the Helmholtz operator. HFH is here applied to flexural waves 
in two dimensions [1591 for which the governing equation is a fourth order equation 



V\ 



Q}u 



0; 



(11.32) 



assuming constant material parameters. Such a thin plate can be subject to point, or line, con- 
straints and these are common place in structural engineering. 

In two dimensions, only a few examples of constrained plates are available in the liter- 
ature: a grillage of line constraints as in Il60l that is effectively two coupled one dimensional 
problems, a periodic line array of point supports ["611 raises the possibility of Rayleigh-Bloch 
modes and for doubly periodic point supports there are exact solutions by [62] (simply sup- 
ported points) and by [63] (clamped points); the simply supported case is accessible via Fourier 
series and we choose this as an illustrative example that is of interest in its own right; it is shown 
in figure ll.ll[ a). In particular the simply supported plate has a zero-frequency stop-band and 
a non-trivial dispersion diagram. It is worth noting that classical homogenization is of no use 
in this setting with a zero frequency stop band. Naturally waves passing through periodically 
constrained plates have many similarities with those of photonics in optics. 

We consider a double periodic array of points at xi =2n\,X2 = 2^2 where w = (with the 
first and second derivatives continuous) and so the elementary cell is one in |.)Ci| < 1, |x2| < 1 
with M = at the origin (see Figure 11.11[ ); Floquet-Bloch conditions are applied at the edges 
of the cell. 
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Applying Bloch's theorem and Fourier series tlie displacement is readily found [[62l as 

/ N / \ v exp(— «;rN-x) ^^^ ^^^ 

m(x) =exp(i»f -x) ) -^ — -,r^ — — =-, (11.33) 

where N = {ni,n2), and enforcing the condition at the origin gives the dispersion relation 

Z)(Ki,K-2,J^) = y 77 -^ z 7^7^ ^ = 0, (11.34) 

In this two dimensional example a Bloch wavenumber vector K = (k"i, K2) is used and 
the dispersion relation can be characterised completely by considering the irreducible Brillouin 



zone YXM shown in figure 11.11 



The dispersion diagram is shown in figure 11.12; The singularities of the summand in 



equation ( 1 1 .34 1 correspond to solutions within the cell satisfying the Bloch conditions at the 
edges, in some cases these singular solutions also satisfy the conditions at the support and are 
therefore true solutions to the problem, a similar situation occurs in the clamped case considered 



using multipoles in [63] . Solid lines in figure 11.12 label curves that are branches of the disper- 



sion relation, notable features are the zero-frequency stop-band and also crossings of branches 
at the edges of the Brillouin zone. Branches of the dispersion relation that touch the edges of 
the Brillouin zone singly fall into two categories, those with multiple modes emerging at a same 
standing wave frequency (such as the lowest branch touching the left handside of the figure at 
M) and those that are completely alone (such as the second lowest branch on the left at M). 

The HFH theory can again be employed to find an effective PDE entirely upon the long- 
scale that describes the behaviour local to the standing wave frequencies and the details are in 



[|59l . the asymptotics from the effective PDE are shown in Fig. |11.12| as the dashed lines. 

11.3 High-contrast homogenization 

Periodic media offer a convenient tool in achieving control of electromagnetic waves, due to 
their relative simplicity from the point of view of the manufacturing process, and due to the 
possibility of using the Floquet-Bloch decomposition for the analysis of the spectrum of the 
wave equation in such media. The latter issue has received a considerable amount of inter- 
est in the mathematical community, in particular from the perspective of the inverse problem: 
how to achieve a given spectrum and/or density of states for the wave operator with periodic 
coefficients by designing an appropriate periodic structure? While the Floquet-Bloch decom- 
position provides a transparent procedure for answering the direct question, it does not yield a 
straightforward way of addressing the inverse question posed above. 

One possibility for circumventing the difficulties associated with the inverse problem is 
by viewing the given periodic structure as a high-contrast one, if this is possible under the val- 
ues of the material parameters used. The idea of considering high-contrast composites within 
the context of homogenization appeared first in the work by Allaire IIT6II . which discussed the 
application of the two-scale convergence technique (Nguetseng [8]) to classical homogeniza- 
tion. A more detailed analysis of high-contrast composites, along with the derivation of an 
explicit formula for the related spectrum, was carried out in a major study by Zhikov [181 . One 
of the obvious advantages in using high-contrast composites, or viewing a given composite as 
a high-contrast one, is in the mere existence of such formula for the spectrum. In the present 
section we focus on the results of the analysis of Zhikov, and on some more recent results for 
one-dimensional, layered, high-contrast periodic structures. 
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In order to get an as short as possible approach to the high-contrast theory, we consider 
the equation of electromagnetic wave propagation in the transverse electric (TE) polarisation, 
when the magnetic field has the form (0,0,//), in the presence of sources with spatial density 

/(x): 

-div(e'^)"^(x/T])V//(x) = (o2//(x)+/(x), x g a C M^ (11.35) 

where we normahse the speed of light c to 1 for simplicity, which amounts to taking £o/io = 1 in 



section 1 1.3 and where the magnetic permeability is assumed to be equal to unity throughout 
the medium {i.e. jj. = jIq), and the function /(x) is assumed to vanish outside some set that 
has positive distance to the boundary of Q.. The inverse dielectric permittivity tensor (e^)^^ (y) 
is assumed in this section, for simplicity, to be a scalar, taking values rj^I and /, respectively, 
on [0, 1] ^-periodic open sets Fq and Fi, such that FqU Fi =M}. Here 7 is a positive exponent 
representing a "contrast" between material properties of the two components of the structure 
that occupy the regions Fq and Fi. In what follows we also assume that Fq fl [0, 1]^ has a finite 
distance to the boundary of the unit cell [0, 1]^, so that the "soft" component Fq consists of 
disjoint "inclusions", spaced [0, 1] ^-periodically from each other, while the "stiff" component Fi 
is a connected subset of E?. The matrix e^ represents the dielectric permittivity of the medium 
at a given point, however the analysis and conclusions of this section are equally applicable to 
acoustic wave propagation, which is the context we borrow the terms "soft" and "stiff" from. 
The assumed relation between the values of dielectric permittivity £^ (in acoustics, between the 
"stiffnesses" ) on the two components of the structure is close to the setting of what has been 
described as "arrow fibres" in the physics literature on electromagnetics, see e.g [[64||. 

A simple dimensional analysis shows that if CO ~ 1 then the soft inclusions are in reso- 
nance with the overall field if and only if 7 = 2, which is the case we focus on henceforth. 



The above equation (11.35) describes the wave profile for a TE-wave in the cylindrical 
domain ^ x M domain, and it is therefore supplied with the Neumann condition dH/dn = 0[j 
on the boundary of the domain and with the Sommerfeld radiation condition dH /d \x\ — i(0H — 
o{\x\^^) as \x\ —7- 00. 

In line with the previous sections, we apply the method of two-scale asymptotic expan- 
sions to the above problem, seeking the solution H = H{x\,X2) = //(x) in the form (see also 
( [mi in Section [TTX2I ) 



//(x)=//o(x,x/t]) + t]//i(x,x/t])-Ft]2//2(x,x/t]) + ..., (11.36) 

where the functions involved are [0, 1] ^-periodic with respect to the "fast" variable y = x/rj. 
Substituting the expansion ( 11.36[ ) into the equation (11.35) and rearranging the terms in the 



resulting expression in such a way that terms with equal powers of 7] are grouped together, we 
obtain a sequence of recurrence relations for the functions //^, fc = 0, 1, ..., from which they are 
obtained sequentially. The first three of these equations can be transformed to the following 
system of equations for the leading-order term //('^)(x,y) = u{x) +v(x,y), x G Q, y G [0, 1]^ : 

-diveholnVi<(x) = (0^(u{x) + f v{x,y)dy) +/(x), xe^, (11.37) 

-Ayv(x,y) = ft)2(^(x)+v(x,y))+/(x), 3; G FoH [0, 1]^, v(x,y)=0, 3; G Fj n [0, 1]^. 

(11.38) 



^Neumann boundary conditions i.e. infinite conducting walls is a good model for metals in microwaves, but 
much less so in the visible range of frequencies wherein absorption by metals need be taken into account. Note 
also that in the TM polarization case, when the electric field takes the form (0,0, £), our analysis applies mutatis 
mutandis by interchanging the roles of e and jj., H and E, and Neumann boundary conditions by Dirichlet ones. 
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These equations are supplemented by the boundary conditions for the function m, of the same 
kind as in the problems with finite r\. For the sake of simplifying the analysis, we assume 
that those inclusions that overlap with the boundary of Q are substituted by the "main", "stiff" 
material, where (e^)^^ = /. 



In the equation (11.37), the matrix Chom is the classical homogenization matrix for the 
perforated medium eF\, see Section above. However, the properties of the system ( 11.37| )- 
( 11.38[ ) are rather different to those for the perforated-medium homogenised limit, described by 
the equation —d.\v£^^^u{x) = tt)^a(x) +/(x). As we shall see next, the two-scale structure of 



( 1 1 .37 )-( 11.38) means that the description of the spectra of the problems ( 1 1 .35 ) in the limit as 
T] — )■ diverges dramatically from the usual moderate-contrast scenario. 

The true value of the above limiting procedure is revealed by the statement of the conver- 
gence, as ?7 — ;■ 0, of the spectra of the original problems to the spectrum of the limit problem 



described above, see [18] and by observing that the spectrum of the system ( 1 1.37 )-( 11-37) is 



evaluated easily as follows. We write an eigenfunction expansion for v(x,y) as a function of 

>'GFon[0,l]2: 



v(x,y) = J^c^(x)v/;t(y), 



(11.39) 



yt=0 



where 1/4 are the (real-valued) eigenfunctions of the Dirichlet problem — Ai//^ = AytV4, j G Fi fl 
[0, 1]^, arranged in the order of increasing eigenvalues Ayt, ^ = 0, 1, ... and orthonormalised ac- 
cording to the conditions /f,n[o.i] 2 IV^/^(y )P^y = 1^ ^ = 0, 1, ..., and /Fonlo,!]^ Wk{y)Wt{y)dy = 0, 
fc 7^ /, fc, / = 0, 1,... Substituting (11.39) into ( |11.38 ), we find the values for the coefficients c^t, 
which yield an explicit expression for v(x, y) in terms of the function m(x) : 

v(x,y)=.(ft)Vx)+/(x))£(/ Wk{y)dy){h-(o^r'Wk{y). 



Finally, using the last expression in the first equation in ( |11.37| ) yields an equation for the 
function u only: 



-diveho!nV"(x) = /3(w') (i<(x) + co-V(x)), x g ^, 
where the function /3, which first appeared in the work [fT8l . is given by 

/3(C02) = C02fl + «2£(^/ VA,(y)Jy)'(A,-«2)-l 



(11.40) 



(11.41) 



The equation (11.40) is supplemented by appropriate boundary conditions and/or condi- 
tions at infinity, which are inherited from the r] -dependent family, i.e. the Neumann condition 
at the boundary points x G <9Q and the radiation condition when |x| — > oo. Clearly, the spec- 
trum of this hmit problem consists of those values of O}^ for which j8 (tt)^) is in the spectrum of 
the operator generated by the differential expression — divgj^J^V subject to the same boundary 
conditions. For example, for the problem in the whole space M^ (describing the behaviour of 
TE-waves in a 3D periodic structure that is invariant in one specified direction) this procedure 
results in a band-gap spectrum shown in Fig. |11.13[ The end points of each pass band are 



found by a simple analysis of the formula ( 1 1 .41 ): the right ends of each pass band are given by 
those eigenvalues X]^ of the Dirichlet Laplacian on the inclusion Fq fl [0, 1]^ that possess at least 
one eigenfunction with non-zero integral over Fq fl [0, 1]^ (otherwise the corresponding term in 
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Figure 11.13: The plot of the function /3 describing the spectrum of the problem {11.37\-{11.38\ subject to the 
boundary conditions. The stop bands for the pro blem i n the whole space K are indicated by the red intervals of 
the horizontal axis. The spectra of the problems \11.35\ considered in the whole space converge, as r\ ^ Q, to the 
closure of the complement of the union of the red intervals in the positive semiaxis. 



( 1 1 .41 1 vanishes), while the left ends of the pass bands are given by solutions to the polynomial 
equation of infinite order /3 (ft)^) = 0. These points have a physical interpretation as eigenvalues 
of the so-called electrostatic problem on the inclusion, see [l23ll . 

As in the case of classical, moderate-contrast, periodic media, the fact of spectral conver- 
gence offers significant computational advantages over tackling the equations ( 11.35[ ) directly: 
as rj — !• the latter becomes increasingly demanding, while the former requires a single numer- 
ical procedure that serves all r] once the homogenised matrix Chom and several eigenvalues X^ 
are calculated. A significant new feature, however, as compared to the classical case, is the fact 
of an infinite set of stop bands opening in the limit as T] — )•, which are easily controlled by the 
explicit description of the band endpoints. This immediately yields a host of applications of 
the above results for the design of band-gap devices with prescribed behaviour in the frequency 
interval of interest. 



The theorem on spectral convergence for problems described by the equation (11.35) is 
proved in [fTSl under the assumption of connectedness of the domain F\ occupied by the "stiff" 
component, via a variant of the extension procedure from Fi to the whole of M^ for function 
sequences whose energy scales as r\^^ (or, equivalently, finite-energy sequences for the oper- 
ator prior to the rescaling x/t] = y). In the more recent works [|24ll . [|25], this assumption is 
dropped in a theorem about spectral convergence for a general class of high-contrast operators. 



via a version of the two-scale asymptotic analysis akin to (11.36), for the Floquet-Bloch com- 
ponents of the resolvent of the original family of operators following the re-scaling x/t] = y. In 
particular, in Il24l a one-dimensional high-contrast model is analysed, which in 3D corresponds 
to a stack of dielectric layers aligned perpendicular to the direction of the magnetic field. Here 
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the procedure described above for the 2D grating fails to yield a satisfactory limit description as 
7] — > 0, i.e. a description where the spectra of problems for finite r] converge to the spectrum of 
the limit problem described by the system ( 11.37[ )-( 11.38) as 7] — ;• 0. A more refined analysis 
of the structure of the related T] -dependent family results in a statement of convergence to the 
set described by the inequalities 



1< 



ha-l5 + l)\^sm(yX{a-l3)]+cos(yX{a-l5)) 



< 1. 



(11.42) 



where a and /3 denote the end-points of the inclusion in the unit cell, i.e. Fq fl [0, 1]^ = (cu,/3) x 
[0,1]. 

Similarly to the spectrum of the 2D high-contrast problem, described by the function /3 , 



the limit spectrum of the ID problem has a band-gap structure, shown in Fig. 11.14 however 
the description of the location of the bands is different in that it is no longer obtained from 
the inequality /3 > 0, where /3 is the ID analogue of (11.41). Importantly, the asymptotic 



behaviour of the density of states function as r/ — ;■ is also very different in the two cases. 



One can show that the family of resolvents for the problems ( 1 1 .35 ) converges, up to a suitable 



unitary transformation, to the resolvent of a certain operator whose spectrum is given exactly 



by (|11.42[), see f25\. The rate of convergence is rigorously shown to be 0(7]), as is anticipated 



by the expansion (11.36) 




Figure 11.14: Tlie square root of t lie limit spectrum for a ID liigli-contrast periodic stac Jc, in T E polarisation. 
Tlie oscillating solid line is the grapli of tlie function f{(0) = cos{(o/2) — a)sin(a)/2)/4 in { 11.42\ with a = 1/4, 
j8 = 3/4. The square root of the spectrum is the union of the intervals indicated by bold lines. 
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The above ID result is generalised to the case of an oblique incidence of an electromag- 
netic wave on the same 3D layered structure. Suppose that X2 is the coordinate across the stack. 
Then, assuming for simplicity that the wave vector (x, 0, 0) is parallel to the direction x\ , it can 
be shown that all three components of the magnetic field are non-vanishing, with the magnetic 
component H — H^ satisfying the equation 

-((£^)-i(x/r])//'w)'=(a)2_(£'7)-i(^/r7)x2)//(x), 
subject to the same boundary condi tions a s before. The modified limit spectrum for this family 



is given by those CO for which (cf. (11 .42 1) 



-l<-ia-l5 + l)(co-—]sm(VX{a-l5))+cos(yX{a-l5))<l, co > 0, (11.43) 
where, as before, a and /3 describe the "soft" inclusion layer in the unit cell, see f24]. The 



set of CO described by the inequalities ( 1 1 .43 1 is similar to that shown in Figure 11.14 the only 



significant difference between the two cases being a low-frequency gap opening near CO = for 



(11.43). 



11.4 Conclusion and further applications of grating theory 

To conclude this chapter, we would like to stress that advances in homogenization theory over 
the past forty years have been fuelled by research in composites [36J . The philosophy of the ne- 
cessity for rigour expressed by Lord Rayleigh in 1892 concerning the Lorentz-Lorenz equations 
(also known as Maxwell- Gamett formulae) can be viewed as the foundation act of homogeniza- 
tion: 'In the application of our results to the electric theory of light we contemplate a medium 
interrupted by spherical, or cylindrical, obstacles, whose inductive capacity is different from 
that of the undisturbed medium. On the other hand, the magnetic constant is supposed to re- 
tain its value unbroken. This being so, the kinetic energy of the electric currents for the same 
total flux is the same as if there were no obstacles, at least if we regard the wavelength as in- 
finitely great.' In this paper, John William Strutt, the third Lord Rayleigh [29], was able to 
solve Laplace's equation in two dimensions for rectangular arrays of cylinders, and in three- 
dimensions for cubic lattices of spheres. The original proof of Lord Rayleigh suffered from 
a conditionally convergent sum in order to compute the dipolar field in the array. Many au- 
thors in the theoretical physics and applied mathematics communities proposed extensions of 
Rayleigh's method to avoid this drawback. Another limit of Rayleigh's algorithm is that it does 
not hold when the volume fraction of inclusions increases. So-called multipole methods have 
been developed in conjunction with lattice sums in order to overcome such obstacles, see e.g. 
[l30l for a comprehensive review of these methods. In parallel to these developments, the quasi- 
static limit for gratings has been the subject of intensive research, one might cite [31] and ll32l 
for important contributions in the 1980s, and [[33ll for a comprehensive review of the modern 
theory of gratings, including a close inspection of homogenization limit. 

Interestingly, in the pure mathematics community, Zhikov's work on high-contrast ho- 
mogenization ifTSl has had important applications in metamaterials, with the interpretation of 
his homogenized equations in terms of effective magnetism first put forward by O'Brien and 
Pendry [[65ll . and then by Bouchitte and Felbacq fi66l . although these authors did not seem to 
be aware at that time of Zhikov's seminal paper [18]. In order to grasp the physical importance 



of ( 1 1.40)-( 1 1.41 ), we consider the case of square inclusions of sidelength a = d/2, where d 
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Figure 11.15: Superlens application of grating: (a) A time liarmonic source at frequency 0.473 display s an image 
througfi a square array of square inclusions; (b) Effective magnetism versus frequency using \11.44 \ for square 
inclusions of relative permittivity 100 witli sidelengtii a — 0.5d in matrix of relative permittivity 1 {grating pitcli 
d = O.lj; Negative values oftlie effective magnetism are in tlie frequency region [0.432,0.534]. 



is the p itch of a bi-periodic grating. The eigenf unctions are \f/, 
(11.41 1 and the corresponding eigenvalues are k^„^ 



n^irp- -\-m^ 



(y) = 2sin(n;r3;i)sin(n7ry2) in 
The right-hand side in the 



homogenized equation ( |11.40[ ) can then be interpreted in terms of effective magnetism: 



l^homK'^) — 1 T~ 



64a^ 



{n,m)odd 



n^m^iklJa^-k^)- 



(11.44) 



This function can be computed numerically for instance with Matlab and demonstrates that 
negative values can be achieved for /i/j^^ near resonances, see Fig. |11.15|^b). This allows for 



superlensing via negative refraction, as shown in Fig. 1 1.15 'a). 

Finally, we would like to point out that high-order homogenization techniques [67] sug- 
gest that most gratings display some artificial magnetism and chirality when the wavelength is 
no longer much larger than the periodicity [|68]| . We hope we have convinced the reader that 
there is a whole new range of physical effects in gratings which classical, high-frequency and 
high-contrast homogenization theories can capture. 
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